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Abstract 

In this paper the concept of wave-guiding medium, previously intro- 
duced for planar structures, is defined for the spherically symmetric case. 
It is shown that a quarter-wavelength layer of such a medium serves as 
a transformer of boundary conditions between two spherical interfaces. 
As an application, the D'B'-boundary condition, requiring vanishing of 
normal derivatives of the normal components of D and B field vectors, is 
realized by transforming the DB-boundary conditions. To test the theory, 
scattering from a spherical DB object covered by a layer of wave-guiding 
material is compared to the corresponding scattering from an ideal D'B' 
sphere, for varying medium parameters of the layer. 

1 Introduction 

Electromagnetic problems are often mathematically defined in terms of differ- 
ential equations and boundary conditions which appear on surfaces limiting the 
region of interest. Here we may separate two different ways to connect mathe- 
matical boundary conditions and a physical structure. In the analytical way one 
starts from a given physical problem involving a structure of electromagnetic 
media and tries to replace existing medium interfaces by boundary conditions 
for the purpose of simplifying the mathematical problem. In such cases the 
boundary conditions are approximative. In the synthetic way a desired struc- 
ture is mathematically designed in terms of boundary conditions, which raises 
the question how to realize them in terms of physical media. This latter way is 
of concern here. 

In [1] it was shown that any given impedance boundary, defined by boundary 
conditions of the type 

u z xE = l s H (1) 
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on a planar surface z — with unit normal u z , can be realized by an interface 
of a wave-guiding medium defined by medium dyadics of the form 



e = e z u z u z + e t , (2) 

p = fi z u z u z +%, (3) 

when the axial parameters grow without limit as e z — > oo, /U 2 — > oo. It was 
shown that the dyadics et and /7 t transverse to the medium axis can be designed 

to realize the surface impedance dyadic Z s . 

As another example, the DB and D'B' boundary conditions on a plane [2, 3], 
respectively defined in terms of normal field components as 

u z • D = 0, u z • B = 0, (4) 

and 

d z (u z -T>) = 0, 8 z (u z -B)=0, (5) 

have been shown to be realizable in terms of physical structures. In fact, already 
in 1959, it was found that the DB boundary can be realized by the interface 
of a uniaxially anisotropic medium (2) and (3) satisfying e z — > and [i z — > 
[4]. The DB medium has recently found application in the design of cloaking 
structures [5, 6, 7, 8]. 

For the planar D'B' boundary, originally introduced in [3], a realization was 
found only recently. In [9] it was shown that a layer of wave-guiding medium, 
called the quarter-wave transformer, upon a planar DB boundary produces D'B' 
boundary conditions at its interface. 

The purpose of the present paper is to study the properties of the wave- 
guiding medium of spherical symmetry in order to extend the previous results 
from planar boundaries to spherical boundaries. The DB and D'B' conditions 
for the spherical boundary take the respective form [3] 



and 



u r • D = D r = 0, 
u r • B = B r = 0, 

V • (u r u r • D) = d r {r 2 D r ) 

V • (u r u r • B) = d r (r 2 B r ) 



(6) 
(7) 

= 0, (8) 
= 0. (9) 



It was shown in [10] that objects with cither DB or D'B' boundary conditions 
with a certain symmetry in their geometry have no backscattering, i.e., they are 
invisible for the monostatic radar. 
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2 Fields in spherical wave-guiding medium 

Let us consider the spherically symmetric anisotropic medium defined by 

e = e r u r u r + e t (\ — u r u r ), (10) 

Jl = /i r U r U r + /j, t (\ - u r u r ), (11) 

where e r ,e t ,/j, r and fi t are constant parameters. The Maxwell equations can be 
represented in spherical coordinates as [11] 

j I \i r rug rsin6>u v \ 

— — det [ d r d e d v = -jiS • H, (12) 

r sm0 \ E r rEg r sin 9E lp J 

^ I u r rug rsm9u ip 

det ( d r dg d v | = juM ■ E. (13) 
H r rHg r sin 9H V 



r 2 sin 9 

Expressing the fields in their radial and transverse components as 



E = —D r + ugEg + UtpEy, (14) 
e r 

H = — B r + UgHg + U v H v , (15) 

and assuming that the radial components of the medium parameters satisfy the 
wave-guiding medium condition 

e r — > oo, ji r — > oo, (16) 

we have E r — > and H r —¥ 0, whence the Maxwell equations for the transverse 
field components are simplified to 

dr{rE v ) = jktVtrHg, (17) 

d r {rE e ) = -jk t VtrH v , (18) 
Vtd r (rH v ) = -jk t rE g , (19) 
Vtd r {rH e ) = jk t rE v , (20) 



where we denote 



UyJ^ttt, Vt = y/vt/et- (21) 



After elimination of fields, the equations become 



(d 2 r + k 2 t )rE v = 0, (22) 



v t J, J-/^ 
\1 i 1.2 \ 



(9 r 2 + ki)rEg = 0, (23) 
(d 2 r + k 2 )rH v = 0, (24) 
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(d 2 + k 2 )rH v = 0, (25) 

It is noteworthy that there is no restriction to the 9 and ip dependence for the 
fields. Thus, the four field components satisfying (22) - (25) can be expressed 
in the form 

F(r) = F + {9,<p)e- jktr + F_(9,tp)e jktr . (26) 

Let us omit the dependence on 9 and ip in the notation. 

After solving the transverse field components, the radial fields B r and D r 
can be obtained from the Maxwell equations as 

B r = ^— (d e (am6E v ) - d v E e ), (27) 
wr sin 9 

D r = 3 — {d e {sm 9H V ) - d^Hg). (28) 

lot sin 

3 Field solutions 

The solutions for the transverse field components can be expressed as 

E v (r) = ^(A+e- jktr + A_e jk * r ), (29) 

Vt H e (r) = --(A+e-^ r - A_e jk * r ), (30) 
r 

E e {r) = ^(C+e- jk * r + C-e jk * r ), (31) 

rhH v (r) = -(C+e- jktr - C_e jk * r ), (32) 
r 

where the coefficients A± , C± are functions of 9 and <p>. The radial fields can be 
obtained by substituting these in (27) and (28). In view of the D'B'-boundary 
conditions (8) and (9), let us define modified radial field quantities as 

B m (r) = r 2 B r (r) 

= -^-Ade{sm9{rE v )) - d^rEg)), (33) 
w sin 9 

D m (r) = r 2 D r (r) 

= J —(de(sm6(rH v ))-d v (rHg)). (34) 

to sin 

The modified fields have the general form 



B m (r) = B + e-J ktr + B_ e i ktr , (35) 
D m (r) = D+e- jk * r + D_e jk * r , (36) 
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with coefficients depending on those of the transverse fields as 

B± = -^—(de(sm0A±) - d v C±), (37) 
u sin 6 

D± = T J —-(do(sm6C ± ) + d v A ± ). (38) 

u>r]t sin 8 

The quantities B mi D m propagate radially for any transverse coordinate 9, ip 
just like waves in a linear transmission line. Forming two new functions involving 
radial derivatives as 

B' m {r) = ld r B m = B+e-^ r - B_e^ r , (39) 

D' m (r) = ld r D m = D + e- jk < r - D_e jk < r , (40) 
the coefficents B±,D± can be solved from (35), (36), (39) and (40) as 

B ± = \{B m {r)±B' m {r))e ± ^, (41) 

D ± = l -(D m (r)±D' m (r))e ± ^\ (42) 

Let us now consider the region between two spherical surfaces, a < r < b. 
Assuming values of B m ,D m ,B' m and D' m known on r = a we can find their 
values on r — b as 

B m (b) \ _ ( B m (a) 



D m (b) \_ K fD T {a) . (44) 



D' m (b) J \ D' m {a) 



through the matrix 



£ = ( cosk t (b-a) -jsmk t (b-a)\ . . 

\ —j sin k t (b — a) cosfc t (6 — a) J 

Similarly, we can find the following relations between the transverse field 
components, 

bE v {b) \_ ( aE v (a) \ 

b m Hg{b) J' \ a Vt Hg(a) )> ^> 

bE e {b) \_ r ( aE e (a) \ 
b Vt H v (b) J ^\ arhH v {a) J ' 



with 



£ = . cosk t (b-a) jsmkt(b-a) . 
' j sin k t (b — a) cos k t (b — a) ' ' 
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4 The quarter-wave transformer 

Let us now assume that the distance between the two spherical surfaces b — a 
satisfies 

kt(b-a) =tt/2, (49) 

whence cos k t (b — a) = and sin k t (b — a) = 1. Defining the wavelength A t by 
fcjAt = 27r, we have b — a = At/4. In this case the layer of wave-guiding medium 
has quarter wavelength thickness and the fields at the two surfaces obey the 
relations 

B m (b) = -jB' m (a), B' m (b) = -jB m (a), (50) 

D m (b) = -jD' m (a), D' m (b) = -jD m (a), (51) 

bE v (b)^ja Vt H e (a), br] t H e (b) = jaE v (a), (52) 

bE e (b) = -jar] t H v (a), br] t H v {b) = -jaE e (a). (53) 

The quarter- wavelength layer serves as a transformer of boundary conditions on 
r = a to other boundary conditions on r = b. Let us list a few examples. 

• DB-boundary to D'B'-boundary 

B m (a)=0 => B' m (b) = 0, (54) 
D m (a)=0 => D^(6) = 0. (55) 

• D'B'-boundary to DB-boundary 

B' m (a)=0 => B m (6) = 0, (56) 
#;»=0 => D m (6) = 0. (57) 

• PEC boundary to PMC boundary 

E v (a)=0 => ifj(6) = 0, (58) 
£ e (a)=0 => i? ¥ ,(6) = 0. (59) 

• PMC boundary to PEC boundary 

H v (a) = => E e (b) = 0, (60) 
#o(a) = => ^(6) = 0. (61) 

• PEMC boundary to PEMC boundary [12] 

^(a) + M^(a)=0 => -1-^(6) =o, (62) 



H e {a)+ME e {a)=0 => H v (6) - -^E v {b) = 0. (63) 

M?7 t 
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The spherical quarter-wave transformer is quite similar to the previously 
studied planar quarter- wave transformer [9]. In particular, the D'B'-boundary 
conditions can be realized on the spherical surface by applying the transformer 
layer upon a DB boundary whose realization is previously known. 

Because the parameters et and m may be freely chosen, the distance between 
the two spherical surfaces can be made as small as we wish by choosing large 
values for e t , /i t . However, we must take care that the parameters e r , fj, r must be 
larger by an order of magnitude, because et/e r and /it//i r must be small. When 
these conditions are met, the transformer layer can be made a thin sheet on the 
sphere of radius a. 

5 Numerical examples 

As an example, let us consider plane- wave scattering from a layered sphere 
simulating the sphere on which a D'B' boundary is forced. The sphere is located 
in free space defined by parameters e ,yu . In the following, we compare the 
scattering behavior of an ideal D'B' sphere and its different material realizations. 

5.1 Volume integral equations 

The scattering fields are found by first solving the equivalent polarization cur- 
rents in the whole sphere (region V) from integral equations. The total time- 
harmonic electric E and magnetic H fields can be expressed via the volume 
equivalence principle [13] as 

E(r) = E mc (r) + — (VV + kf\) • S{3) - V x 5(M), (64) 

H(r) = H mc (r) + V x 5(J) + — (VV + kf\) • «S(M), (65) 

where E mc and H mc are the incident electric and magnetic fields, and k = 
ojy/e [i . The volume integral operator in (64) and (65) is defined by 

5(F)- J G (r,r')F(r')oV', (66) 
v 

where G is the free-space Green function. The equivalent electric and magnetic 
polarization currents are defined as 

J(r) = jw(f(r) - ej) ■ E(r) 

M(r) = jw(p(r) - MoT) • H(r). 

By applying the volume equivalence principle, definitions of the volume 
equivalent currents J and M, vector identities, and the fact that the operator 
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(66) satisfies the Helmholtz equation, we obtain the following volume integral 
equation formulation, 

-jcje f E • E mc (r) = -1 • J + W E ■ (V x V x 5(J)) - jwe W E • V x 5(M), (68) 

-j^ f M -H mc (r) = -p-J + jcj Mo f M -Vx5(J)+f M -(Vx Vx5(M)), (69) 
with the contrast dyadics defined by 

e TE = e - e T, HqTm = ~p- MoT- (70) 

The formulation (68), (69) can be discretized using Galerkin's method with 
piecewise constant vector basis and testing functions. The hypersingularity of 
the kernel is reduced by moving one derivative into the testing function through 
integrating by parts. The remaining derivatives are then moved to operate on 
the Green function. 

5.2 Scattering from a layered sphere 

Let us now consider scattering from a layered sphere of radius b = 1 m, for an 
incident plane wave with free-space wavelength A Q = 3.2 m, propagating in the 
direction of positive z axis (9 = 0). Fig. 1 presents two possible realizations 
for the D'B' sphere with transverse parameters either \i t — Ho, c t = e or 
fit = 2/i Q , e t = 2e D . In order to obtain DB conditions at the surface r = a, 
the permittivity and permeability of the inner sphere are set to zero [14]. In 
Fig. 1 (left) the thickness of the quarter-wave transformer layer is b — a = 0.8 m 
since X t = \ = 3.2m. In Fig. 1 (right), the reduced wavelength becomes 
A* = A G /2 = 1.6m and therefore the layer thickness equals b — a = 0.4m. 

Fig. 2 shows the scattering cross section (SCS) of the coated sphere in Fig. 1 
(left), as calculated by the volume integral equation method. The calculated 
SCS of the sphere of Fig. 1 (right) is presented in Fig. 3. We have used dif- 
ferent values of radial permittivity and permeability in order to get an idea on 
how large radial components of permittivity and permeability are required in 
approximating the D'B' boundary. Note that, due to the self-dual character 
of the structure, the scattering pattern is rotationally symmetric, i.e. E and H 
planes have the same patterns. 

Clearly, radial parameter values e r = 10e o , \i r — 10^ o yield a poor approxi- 
mation to the quarter- wave transformer, whence the sphere markedly deviates 
from one with the D'B' boundary. However, for larger values of radial permit- 
tivity and permeability, the scattering cross sections of the two coated spheres 
approach to that of the ideal D'B' sphere. The ideal D'B' scattering is com- 
puted from an exact Mie scattering solution [10]. From the numerical results it 
appears that, for larger values for the transverse parameters et,Ht, lower values 
of the radial parameters e r , \i r can be accepted for the realization. Furthermore, 
one may note that, for all parameter values used in the computation, the layered 
sphere has very small scattering in the backward direction (9 — 180°). This is 
in accord with the theory stating that both DB and D'B' spheres are invisible 
for the monostatic radar [10]. 
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Figure 1: D'B'-boundary conditions on a spherical surface of radius b = 1 m can 
be realized by a quarter-wave transformer layer upon a sphere of radius a with 
DB-boundary conditions. The thickness b — a of the layer can be controlled 
by the components of permittivity and permeability transverse to the radial 
direction in the layer, et>Mt- The two cases are considered here as examples. 

6 Conclusion 

An anisotropic medium with infinitely large axial permittivity and permeabil- 
ity has been called waveguiding medium and its properties in transforming 
impedance-boundary conditions between two planar interfaces has been studied 
in [1]. In [9] it was shown that a layer of waveguiding medium can be also 
used to transform between planar DB and D'B' boundary condition which in- 
volve vanishing of normal components of the D and B fields or their normal 
derivatives. In the present paper the spherical waveguiding medium is defined 
and its transforming properties are studied. In particular, it is shown that a 
quarter-wavelength layer of waveguiding medium can be used to transform a DB 
boudary to a D'B' boundary. Since it is known that a spherical DB boundary 
can be realized by a medium with vanishing permittivity and permeability [14], 
this gives a means to realize a spherical D'B' boundary which until now has had 
no realization whatever. The theory is verified numerically through volume- 
integral equation approach by considering scattering from a layered spherical 
object and comparing with that from an ideal D'B' sphere. It is seen that the 
realization approaches the ideal case when the radial parameters of the quarter- 
wavelength layer grow large. Such an object is of interest because it has zero 
backscattering like the corresponding DB sphere [10]. 
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